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ABSTRACT 

2-Amino-3-methylimidazo[4,5-f|quinolone (IQ), a 
heterocyclic amine found in coolced meats, 
undergoes bioactivation to a nitrenium ion, whicli 
all<ylates guanines at botli the C8-dG and A/^-dG 
positions. The conformation of a site-specific 
A/^-dG-IQ adduct in an oligodeoxynucleotide 
duplex containing the iterated GG repeat restriction 
site of the Nar\ endonuclease has been determined. 
The IQ moiety intercalates, with the IQ H4a and GH3 
protons facing the minor groove, and the IQ H7a, 
H8a and H9a protons facing the major groove. The 
adducted dG maintains the anti-conformation about 
the glycosyl bond. The complementary dC is 
extruded into the major groove. The duplex main- 
tains its thermal stability, which is attributed to 
stacking between the IQ moiety and the 5 - and 
3'-neighboring base pairs. This conformation is 
compared to that of the C8-dG-IQ adduct in the 
same sequence, which also formed a 'base- 
displaced intercalated' conformation. However, the 
C8-dG-IQ adopted the syn conformation placing 
the Watson-Crick edge of the modified dG into 
the major groove. In addition, the C8-dG-IQ adduct 
was oriented with the IQ CH3 group and H4a and 
H5a facing the major groove. These differences 
may lead to differential processing during DNA 
repair and replication. 

INTRODUCTION 

Although browning meats during cooking imparts flavor, 
it also leads to the formation of heterocyclic amines 
(HCAs) such as 2-amino-3-niethylimidazo[4,5-/]quinoline 
(IQ) (1-5). IQ has been produced in cooked meats at 
ppb levels (6,7) and has also been detected in tobacco 



smoke (8). HCAs and their metabolites have been 
isolated from human urine (9). Human exposures to 
HCAs, estimated to be ~60ng/day (10), are modest, but 
are hkely to be involved in cancer etiology (11,12). 

IQ is a potent mutagen (13). It is less prevalent than 
2-amino-l-methyl-6-phenylimidazo[4,5-b]pyridine (PhIP) 
(14), but it is 200-fold more mutagenic in Salmonella re- 
version (Ames) assays (3) and it is an order of magnitude 
more mutagenic than aflatoxin Bi. In these assays (15-18), 
HCAs such as IQ are active in point and frameshift tester 
strains (19). In bacteria, mutations occur primarily at G:C 
base pairs (20,21). IQ is a potent inducer of two-base 
frameshifts in CG repeats. Similar levels of mutations 
are seen in mammahan hprt (22) and ef-2 (23) gene 
assays. Base-pair substitutions are the predominant muta- 
tions observed in mammahan cells (24-26). Sister chroma- 
tid exchange has been observed in rodent cells (27-29). IQ 
induces tumors in the organs of rodents and in the livers of 
monkeys (30-33). Liver, forestomach and lung tumors 
have been observed in IQ treated mice (34), whereas 
liver, intestine, zymbal gland, clitoral gland, skin (35), 
mammary glands, liver and ear ducts tumors have been 
observed in exposed rats (36). The TD50 values for IQ are 
0.7 and 14.7nig/kg/day in rats and mice, respectively (37). 
Human exposures to HCAs have been associated with 
pancreatic (38), colon (39), prostate (40) and breast 
cancers (41,42). 

The genotoxicity of IQ derives primarily from its oxi- 
dation by CYP P450 1A2 to an jV-hydroxylamine (43^7) 
although extra-hepatic CYP P450s oxidize HCAs with 
lower efficiencies (Scheme 1) (48). The A^-hydroxylaniine 
is acetylated by A^-acetyl transferases, particularly NAT2 
(49-51). In humans, the NAT2 fast acetylator polymorph- 
ism correlates with increased genotoxicity and cancer 
(52-54). The nitrenium ion resulting from solvolysis of 
A^-acetoxy-IQ is the ultimate electrophile (26,48). It 
reacts predominately at the C8 atom of guanine, while a 
minor alkylation product is formed at the A^^ atom of 
guanine (55-57). In addition, IQ can be converted to a 
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Scheme 1. The bioactivation of IQ. Cytochrome P450-mediated W-hy- 
droxylation followed by W-acetylation of the hydroxylamine and 
deacetylation forms an electrophilic aryl nitrenium ion. The aryl 
nitrenium ion alkylates guanine in DNA either via the IQ amine 
nitrogen to form the C8-dG-IQ adduct, or alternatively, alkylates 
DNA via the C5 position of the IQ ring to form the W^-dG-IQ adduct. 



reactive and genotoxic A^-nitrosamine, which shows 
similar regioselectivity for DNA alkylatioii (58,59). 

Levels of C8 and A^^-dG-lQ adducts have been 
measured in rat and primate tissues using mass spectrom- 
etry (60,61) and "^^P post-labeling methodology. Turesky 
et al. (62) have monitored C8-dG-IQ adduct formation in 
human hepatocytes using tandem Hquid chromatography- 
electrospray ionization mass spectrometry. Levels range 
from 7 to 26 adducts per 10^ bases. While less abundant, 
the jV^-dG-IQ adduct is more persistent than is the C8-dG- 
IQ adduct, suggesting that it is repaired less efficiently 
(63). The A^^-dG-IQ adduct may therefore play a signifi- 
cant role in the genotoxicity of IQ. 

We have synthesized phosphoramidite reagents of the 
C8- and jV^-dG-IQ adducts in which the Buchwald- 
Hartwig palladium-catalyzed A^-arylation was the key 
C-N bonding-forming step (64-67). These adducts have 
been site-specifically incorporated into the Narl restriction 
sequence, 5'-d(CG'G^CXXC)-3', using automated solid- 
phase synthesis. The G^^ position is a hot spot for two-base 
frameshift deletions in bacterial mutagenesis assays, while 
the G' position is not (68-71). In addition, human DNA 
polymerase (hpol) ri produces two-base deletions when 
replicating past the A^'^-dG-IQ adduct at position G^^, 
in vitro (72). Thus, this sequence provides a platform for 
investigating sequence-specific conformational perturb- 
ation of DNA structure by IQ adducts, in relationship 
to their biological processing (68-71,73,74). Previously, 
we determined the conformation of the C8-dG-IQ 



adduct at the G^ position of the Narl sequence, which 
exhibited a base-displaced intercalated conformation (75). 

Presently, we have determined the conformation of 
the 7V^-dG-IQ adduct in 5'-d(C'T^C^G'^G^C*^X''C**C'' 
A"'T"c'V:5'-d(G'^A'4T'^G'«G''c'**G"'c2°C-^2^A2^ 
G-'*)-3'; X = A^^-dG-IQ adduct (Chart 1). This duplex 
contains the recognition sequence of the Narl restriction 
endonuclease, in which the G^ nucleotide (X"* in the Narl 
sequence and X^ in this study) represents a hot spot for 
two-base deletions. The IQ moiety intercalates, with the 
IQ H4a and CH3 protons facing the minor groove, and the 
IQ H7a, H8a and H9a protons facing the major groove. 
The adducted nucleotide maintains the anti-conformation 
about the glycosyl bond. The complementary dC is 
extruded into the major groove. Nevertheless, the duplex 
maintains its thermal stabihty. This is attributed, in part, 
to stacking between the IQ moiety and the 5'- and 
3'-neighboring base pairs. The base-displaced intercalated 
conformation of the A^^-dG-IQ adduct differs from that of 
the C8-dG-IQ adduct, and may also be compared with 
that of the A^^-dG AAF adduct (76), providing insight as 
to the persistence of the iV^-dG-IQ adduct (63) and its 
processing during replication and repair events. 



MATERIALS AND METHODS 

Sample preparation 

The TV^-dG-IQ-adducted oligodeoxynucleotide 5'-d(CTC 
GGCXCCATC)-3' was synthesized as described (67). 
The complement strand 5'-d(GATGGCGCCGAG)-3' 
was synthesized by Midland Certified Reagents Co. 
(Midland, TX, USA) and purified by anion exchange 
chromatography. HPLC chromatographic utilized a 
Supelcosil LS-18-DB analytical base-deactivated C-18 
column (Sigma-Aldrich Inc., St. Louis, MO, USA), 
using a gradient from 5% to 12% acetonitrile in 
ammonium formate buffer (pH 7), over 35min, or a 
gradient from 5% to 30% acetonitrile in NaH2P04 
buffer (pH 5), over 25min. The oligodeoxynucleotides 
were characterized by negative mode MALDI-TOF spec- 
trometry in a hydroxypicolinic acid matrix. The 
ohgodeoxynucleotides were annealed at 1:1 ratio at 
room temperature in 180 |il buffer containing lOmM 
NaH2P04, 100 niM NaCl and 5nM NasEDTA (pH 7). 

Thermal melting experiments 

UV melting temperatures were collected on Cary 100 Bio 
UV spectrometer using 0.5 OD of duplex in 1ml of 
solution containing 0.1 M NaCl, lOniM NaH2P04 and 
0.05 mM Na2EDTA (pH 7.0). The temperature was 
increased from 25 to 75° C at a rate of 1°C per min. 

NMR spectroscopy 

The A^^-dG-IQ modified and the unmodified duplexes 
were prepared at concentrations of 570 and 810 ^iM, re- 
spectively, and placed into 3 mm diameter micro NMR 
sample tubes (77). The samples were prepared in 0.1 M 
NaCl, 50|.iM NajEDTA and 10 mM NaH2P04 (pH 7.0). 
To observe non-exchangeable protons, the samples were 
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2-Amino-3-methylimida20[4,5-/|quinoline (IQ) N2(2-amino-3-methyl-3H-imidazo[4,5-f|quinol-5-yl)-dG 




C8-[(3-methyl-3H-imidazo[4,5-f|quinolin-2-yl)amino]-dG 3-(deoxyguanosin-N2-yl)-2-acetylaminofluorene Adduct 
E 

5'-Cl t2 c3 g4 g5 c6 c8 c9 aIOtIIcIS-S ■ 

3 ■ _g24a23g22c21c20q19c18g17g16t15a14g13_5 . 

Narl Oligodeoxynucleotide Sequence 

Chart 1. (A) Structure of IQ. (B) Structure of the W^-dG-IQ adduct, showing the numbering of guanine base and IQ protons. (C) Structure of the 
C8-dG-IQ adduct. (D) Structure of the W^-dG-AAF adduct. (E) The duplex containing the Narl sequence, showing the numbering of the nucleotides. 
The W'-dG-IQ adduct is positioned at X^, which corresponds to the G^ frameshift-prone position of the Narl sequence. 



exchanged with D2O. NMR spectra were recorded at 
600 or 800 MHz. The spectra were collected at 15°C; 
NOESY experiments were conducted at mixing times of 
150, 200 and 250ms with a relaxation delay of 1.8s. 
Additional experiments were condticted with a longer re- 
laxation delay to evaltiate NOE distances arising from the 
adenine H2 protons, which typically exhibit longer Ti re- 
laxation values. The data were collected with 512 points 
in the /I dimension and 2048 points in the /2 dimension. 
Chemical shifts were referenced to water. For the 
observation of exchangeable protons, the samples 
were dissolved in 9:1 H20:D20. 'H NMR spectra were 
recorded at 600 or 800 MHz at 15°C. The data were 
collected with 512 points in the tl dimension and 2048 
points in the t2 dimension. A mixing time of 250 ms 
was used. Water suppression was performed using the 
WATERGATE pulse sequence (78). The spectra were 
processed using the TOPSPIN software (Bruker Biospin 
Inc., Billerica, MA, USA). 

NMR experimental restraints 

The spectral data were evaluated using the program 
SPARKY (79). The intensities of NOE cross peaks were 
measured by volume integrations. The bounds for 
overlapped peaks were optimized manually. Noise was 
assigned half the intensity of the weakest peak, and 
motion was assumed to be isotropic. Experimental 
intensities were combined with intensities obtained from 
complete relaxation matrix analysis (CORMA) of starting 



model to generate a hybrid intensity matrix (80,81). The 
intensities were converted to distances with the program 
MARDIGRAS (82), which refined the hybrid intensity 
matrix. Calculations were performed using 150, 200 and 
250 nis mixing time data and 2, 3 and 4 ns isotropic cor- 
relation times. Evaluation of the resulting distance data 
allowed creation of upper and lower bound distance re- 
straints that were used in restrained molecular dynamics 
(rMD) calculations. 

Restrained molecular dynamics calculations 

An unmodified B-DNA model (83) was used as a starting 
structure. The guanine at position was replaced by the 
7V^-dG-IQ adduct using the program INSIGHT II 
(Accelrys Inc., San Diego, CA, USA). Partial charges 
for the iV^-dG-IQ adduct were calculated with the 
B3LYP/6-31G* basis set in GAUSSIAN (84). The 
starting structure was energy minimized for 1000 cycles. 
Simulated annealing protocols (85) used for the rMD cal- 
culations were conducted with the parm99 force field (86), 
using the program AMBER (87). Force constants of 
32 kcal/mol/A^ were apphed for distance and torsion 
angle restraints. The generalized Born model (88) was 
used for solvation. The salt concentration was 0.1 M. 
The molecule was coupled to the bath to control the tem- 
perature during simulated annealing (89). First, calcula- 
tions were performed for 20 ps (20000 steps) and 
recording data every ps by the following protocol: 
during steps 0-1000, the system was heated from 0 to 
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600 K with a coupling of 0.5 ps. During steps 1001-2000, 
the system was kept at 600 K. The system was then cooled 
from 600 to 100 K during steps 2001-18 000 with a 
coupHng of 4ps. Further cooling from 100 to OK 
occurred during steps 18 001-20000 with a couphng of 
1 ps. After initial cycles of refinement a longer 100 ps 
(100 000 steps) calculation was performed by the following 
protocol: During steps 0-5000, the system was heated 
from 0 to 600 K with a coupling of 0.5 ps. During steps 
5001-10000, the system was kept at 600 K. The system 
was cooled from 600 to 100 K during steps 10001-90000 
with a coupling of 4ps. Additional cooUng from 100 to 
OK occurred during steps 90001-100000 with a coupling 
of 1 ps. Structure coordinates were saved after each cycle 
and were subjected to potential energy minimization. 
CORMA (80,81) was used to compare intensities 
calculated from these emergent structures with the 
distance restraints. Helicoidal analysis was performed 
using the CURVES+ web server (90,91). 

RESULTS 

Oligodeoxynucleotide containing the A'^-dG-IQ adduct 

The iV^-dG-IQ adduct was incorporated into 5'-d(CTCGG 
CXCCATC)-3' using automated solid-phase synthesis (67). 
The position of the A^^-dG-lQ adduct was located at 
position X^, corresponding to position G"* in the Narl 
sequence. The modified oligodeoxynucleotide was purified 
by CI 8 reverse phase HPLC and characterized by MALDI- 
TOF mass spectrometry in negative ion mode [m/z 3777.7, 
calcd for (M - H), 3776.6]. Thermal melting (Tn,) profiles of 
0.5 Ajeo units of the IQ-modified duplex were monitored at 
100 niM NaCl (1 ml volume) as a function of temperature 
by absorbance at 260 nm. An unmodified duplex was 
evaluated under the same conditions to provide a basis of 
comparison. The of the modified duplex was 63° C, 
within experimental error of the unmodified duplex 
(Supplementary Figure SI). Thus, the A^^-dG-IQ adduct 
did not reduce the stabiHty of this oligodeoxynucleotide. 
This result differed from our previous report, which had 
indicated that this adduct destabilized this duplex (67). 
Subsequent analysis of the previous sample by mass spec- 
trometry revealed that the complement strand was not 
correct, accounting for the discrepancy. Table 1 lists the 
correct values of the A'^^-dG-IQ adduct at the three pos- 
itions of the Narl sequence. 

NMR 

The modified duplex yielded well-resolved NMR spectra 
with narrow line shapes for the non-exchangeable protons 
at 15°C. The best spectral quality for the exchangeable 
protons was obtained at 5°C. 

Non-exchangeable DNA protons 

The base aromatic and deoxyribose anomeric protons were 
assigned using established procedures (Figure 1) (92,93). 
The intensity of the H8 to X^ HI' NOE was not 
changed in the presence of the adduct, indicating minimal 
change in the conformation of the glycosyl torsion angle. In 
the complementary strand, the intensity of the NOE 



Table 1. Thermal melting temperatures (T^ measurements) of Narl 
duplexes containing the W^-dG-IQ adducts 



Narl N^-dG-lQ modified duplex 


Tm (°C) 


Ar„* (°c) 


5'-CTCXGCGCCATC-3' 


62 


-1 


3'-GAGCCGCGGTAG-5' 






5'-CTCGXCGCCATC-3' 


64 


+1 


3'-GAGCCGCGGTAG-5' 






5'-CTCGGCXCCATC-3' 


63 


0 



3'-GAGCCGCGGTAG-5' 



between C'** HI' and G''' H8 was weakened. The ^^^-dG- 
IQ adduct did not induce breaks in the sequential pattern of 
NOEs between the aromatic base protons and the anomeric 
protons. With the exception of the adduct site, the 
internucleotide NOEs were characteristic of a B-type 
duplex. The adenine H2 protons were assigned based 
upon NOEs to the thymine imino protons of the respective 
A:T base pairs. With the deoxyribose HI' assignments in 
hand, the remainder of the deoxyribose protons was 
assigned from a combination of NOESY and COSY 
data. The assignments of the non-exchangeable DNA 
protons are summarized in Supplementary Table SI. 

Exchangeable DNA protons 

The imino and amino proton regions of the NOESY 
spectrum are shown in Figure 2. The assignments were 
made using established methods (94). The A^^-dG-IQ 
adduct perturbed Watson-Crick hydrogen bonding. At 
the X^:C'^ base pair, the X^ imino proton resonance was 
broadened, probably due to an enhanced rate of exchange 
with water. The amino protons for C'** were not detected. 
No NOE was observed between the X'' and G'^ imino 
protons, perhaps due to the broadening of the X' imino 
proton. The chemical shifts of the X and G'^ imino 
protons were almost isochronous. It was not possible to 
determine if a NOE between these two protons existed. All 
other base pairs were assigned, with the exception of the 
two terminal base pairs C'lG^"* and C'^:G'^. The imino 
protons from the terminal base pairs were exchange 
broadened. Overall, the data suggested that the duplex 
maintained Watson-Crick hydrogen bonding, with the ex- 
ception of the modified base pair (Figure 2B). The assign- 
ments of the exchangeable protons are summarized in 
Supplementary Table S2. 

IQ protons 

The IQ protons, consisting of the CH3 group, the H4a 
proton, and the H7a, H8a and H9a spin system, were 
assigned from a combination of COSY and NOESY 
data (Figure 3). The CH3 resonance was observed at 
3.57 ppm. It displayed an intense NOE to the H4a 
proton, whose resonance was observed at 8.55 ppm. A 
coupling between the H8a proton (8 6.55 ppm) and the 
H9a proton (§ 7.65 ppm) was observed in the COSY 
spectrum. The H8a proton also exhibited an NOE to the 
H7a proton (5 7.6 ppm). The coupling between H8 and 
H7a exhibited weak intensity in the COSY spectrum. This 
was attributed to presence of the nitrogen atom in the 
ring, which broadened the H7a resonance. This effect 
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Figure 1. Expanded plot of the 250 ms NOESY spectrum showing NOEs between the base aromatic and deoxyribose anomeric protons of the N - 
dG-IQ modified duplex. (A) Bases C' to C'^ of the modified strand. (B) Bases G'"" to G^"* of the complementary strand. The spectrum was acquired 
at 800 MHz at 15°C. 




Figure 2. Expanded plots of the NOESY spectrum, showing the NOEs between the exchangeable imino and amino protons of the W'^-dG-IQ- 
modified duplex. (A) Sequential NOE connectivity for the imino protons of base pairs T^:A-^ to T":A'''. (B) NOE connectivity within Watson-Crick 
base pairs and between the imino protons and the amino protons. The lettered cross-peaks are assig ned as follows: a', N^U -^C^ A^Hb; a, 
N^n ^C*" A^Ha; b', G" NIH^C* Ar*Ha; b, G" NIH^C* A/^Hb; c', G" NIH^C** A^Hb; c, G" NIH^C*" A/^Ha; d', G"" NIH^C'' A^^Hb; d, G"* 
NIH^C'' A^Ha; e', G^^ NIH^C' A^Hb; e, G^^ NIH^C' A^Ha; f , G' N1H^C^° A^Hb; f, G^ NIH^C^" A^Ha; g', G" NIH^C^' A^Hb; g, G" 
NIH^C^' A^Ha; h, T'^ NaH^A'" H2; i, T" N3H^A''' H2; j, T^ N3H^A^' H2; 1, X^ W^H^IQ H4a; 2, IQ H7a^X' A^^H; 3, X'** NIH^IQ 
H4a. The spectrum was collected at 800 MHz at 5°C. 



was also observed for the C8-dG-IQ adduct, for which the 
COSY cross-peak between H7a and H8a was only 
observed between 25 and 35°C (95). The IQ amine 
proton was not assigned. 

Chemical shift pevturhations 

The A^^-dG-IQ adduct resulted in localized chemical shift 
perturbations, involving the modified base pair X':C'** 
and the neighboring C'':G'' and C**:G'^ base pairs 



(Figure 4). At the modified X^:C"^ base pair, the H8 
resonance shifted 0.4 ppm downfield relative to the G^ 
H8 resonance in the unmodified duplex. In contrast, the 
C'** H6 and C'** HI' resonances shifted 1 and 0.8 ppm 
downfield, respectively. At the 5'-neighbor C^:G'' base 
pair, the H6 resonance shifted upheld by 0.2 ppm, 
whereas the HI' resonance shifted downfield by 
0.4 ppm. The G" H8 and HI' resonances each shifted 
upheld by 0.4ppm. At the 3'-neighbor C**:G'^ base pair, 
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Chemical Shift Perturbations 



ppm 



Figure 3. Expanded tile plot of the 250 ms NOESY spectrum showing 
the assignments of the IQ ring protons. (A) The IQ H4a proton was 
observed at 8.55 ppm. (B) The IQ H7a proton was observed at 
7.65 ppm. (C) The IQ H9a proton was observed at 7.6 ppm. (D) The 
IQ H8a proton was observed at 6.55 ppm. (E) The IQ CH3 protons 
were observed at 3.57 ppm. The spectrum was collected at 800 MHz at 
15°C. 



the C** H6 resonance shifted downfield by 0.3 ppm, 
whereas the C*^ HI' resonance shifted upheld by 0.2 ppm. 
The G'^ H8 resonance sliifted downfield by 0.1 ppm. The 
G'^ HI' resonance showed negligible chemical shift per- 
turbation. The resonances for the remaining base pairs in 
the duplex also showed neghgible chemical shift perturb- 
ations. In the imino proton region of the spectrum, the 
and G'^ NIH imino resonances, at 11.57 and 11.59 ppm, 
respectively, exhibited upheld chemical shifts of > 1 ppm 
from those of the unmodified duplex, at 13.24 and 
13.16ppm, respectively. The A^^H amine resonance 
was observed at 9.5 ppm. 



NOEs between IQ and DNA 

The CH3, H4a, H7a and H8a protons of IQ exhibited 
NOEs to the C^ G", C'** and G" bases (Table 2). The 
pattern of NOEs involving H9a was difficult to establish 
due to resonance overlap with G'^. The CH3 group 
showed medium strength NOEs to X^ HI' and C*^ H6 
and weak NOEs to HI' and C** H5. It also showed 
an NOE to the G'^ NIH imino proton. The H4a proton 
showed a strong NOE to X'' HI', and medium NOEs to C** 
H5, the X'^ N2H amine proton and the G" NIH imino 
proton and weak NOEs to X'' H2' and X' H2" The H7a 
proton exhibited weak NOEs to G'^ HI', and the X^ N2H 
amine proton. Three NOEs from X^ to C'** were observed; 
these were of medium strength NOE between H9a and C'** 
HI', and H8a and medium strength between C'** H2' and 
H2". The H8a proton showed a medium strength NOE to 
G'^ H3' and weak NOEs to G'^ HI' and G'' H8. Some 
30% of the NOEs from the IQ ring were to protons in the 
complementary strand, whereas another 46% were to 
other IQ protons and protons of the modified base. The 
remaining 24% of the NOEs were to neighbor bases in the 
modified strand. 
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Figure 4. Chemical shift perturbations of the deoxyribose HI' protons 
(gray) and the pyrimidine H6 or purine H8 aromatic protons (black), 
for the 7V^-dG-IQ modified duplex. (A) Nucleotides C'-C'^ in the 
modified strand. (B) Nucleotides G'^^-G^'' in the complementary 
strand. The A5 (ppm) values were calculated as S^ojiUed dupiex^Sunmodi- 
fied duplex. Positive ASvalues represent upheld chemical shift perturb- 
ations. Negative ASvalues represent downfield chemical shift 
perturbations. 



Conformational Refinement 

After the unmodified duplex was constructed using 
B-DNA coordinates (83), the guanine at position G^ was 
replaced by the A^"-dG-IQ adduct. The partial charges for 
the A^'^-dG-IQ adduct are provided in Supplementary 
Figure S2. Potential energy minimization provided an 
energy minimized starting duplex. A total of 329 
distance restraints consisting of 127 inter- and 202 intra- 
nucleotide distances (Table 3) were obtained using the 
program MARDIGRAS (81,82), from 15°C NOESY 
data. Similar distance restraints were obtained if the 
data were collected at 150, 200 or 250 ms mixing times. 
These restraints included 16 DNA-IQ distances. A total of 
49 Watson-Crick hydrogen-bonding restraints were 
applied for all of the base pairs except for the modified 
X'^iC'** base pair. An additional 100 phosphodiester 
backbone and 20 deoxyribose pseudorotation restraints 
for base pairs not proximal to the site of modification 
were obtained from canonical values derived from 
B-DNA (83), consistent with the spectroscopic data 
indicating that the duplex maintained a B-DNA like struc- 
ture. A series of rMD calculations were performed using a 
simulated annealing protocol in which the generalized 
Born solvation model (88) was used, with a salt concen- 
tration of 0.1 M. The emergent structures were subjected 
to potential energy miniinization before further analysis, 
which involved a 100 ps rMD calculation using the 
protocol described above, again followed by potential 
energy minimization. 

The pairwise rmsd analysis of structures emergent from 
the rMD calculations was used to measure the precision of 
the structural refinement. Ten structures were chosen 
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Table 2. Summary of NOEs observed between A''^-dG-IQ(X7) adduct protons and oligodeoxynucleotide protons and their intensities 



IQ proton 


NOEs to oligodeoxynucleotide protons 




CH3 


HI': medium; C** H6: medium; C** HI 


': weak; C** H5: weak 


H4a 


HI': strong; C** H5: medium; X' H2': 


weak; X' H2": weak; X^ N2H: medium; G" NIH: medium 


H7a 


G'' HI': weak; X^ H2; weak 




H8a 


G" H3': medium; C"* H2': medium; C'** 


H2": medium; G'^ HI'; weak; G" H8: weak; G'' H2"; weak 


H9a 


C"* HI': medium 





Table 3. NMR restraints used for the W^-dG-IQ stri 
tions and refinement statistics 



NOE restraints 
Internucleotide 
Intranucleotide 
Total 

Backbone torsion angle restraints 
H-bonding restraints 
Deoxyribose restraints 
Total number of restraints 
Refinement statistics 
Number of distance restraint violations 
Number of torsion restraint violations 
Total distance penalty/maximum penalty (kcal/mol) 
Total torsion penalty/maximum penalty (kcal/mol) 
r.m.s. distances (A) 
r.m.s. angles (°) 

Distance restraint force field (kcal/mol/A^) 
Torsion restraint force field (kcal/mol/deg^) 



;ture calcula- 



127 

202 

329 

100 

49 

20 

498 

56 
50 

2.3/0.187 

2.8/0.177 

0.012 

2.5 

32 

32 




based on the lowest deviations from the experimental 
distance and dihedral restraints (Figure 5). These ex- 
hibited an rmsd of 0.012 A in distances and 2.5° in 
torsion angles (Table 3). There were 56 distance violations 
with a maximum penalty of 0T87 kcal/mol and a total 
distance penalty of 2.3 kcal/mol. There were 50 torsion 
angle violations with a maximum penalty of 0.177 kcal/ 
mol and a total torsion angle penalty of 2.8 kcal/mol. 
The maximum pairwise rmsd distances were 1.1 2 A. 
These structures were averaged and subjected to potential 
energy minimization. 

The accuracy of the refined structures was assessed by 
complete relaxation matrix analyses (80,81), which 
compared intensities calculated from the refined structures 
with the distance restraints (Figure 6). The sixth root 
residual i?'x value of the average structure was 8.4%, 
and the individual values for intra-nucleotide restraints 
(8.5%) and inter-iiucleotide restraints (8.3%) were of 
siinilar magnitudes. This indicated agreement with the 
NOE data. Nucleotide G'^ exhibited a greater i?'x value 
of 17.1%, suggesting that it was not as well-refined. This 
was attributed to several NOEs involving G'^ being 
overlapped with other resonances. The structural statistics 
are summarized in Table 4. 

Conformation of the A'^-dG-IQ Adduct 

The modified nucleotide (X^) remained in the anti-con- 
formation about the glycosyl bond. It was displaced 
toward the major groove. The IQ ring was intercalated 




Figure 5. Superposition of ten potential energy minimized structures 
emergent from the rMD calculations of the A'~-dG-IQ modified 
duplex, using distance restraints from the 250 ms NOESY data. The 
positions of the modified X' nucleotide and the C"* nucleotide in the 
complementary strand are as indicated. The maximum pairwise rmsd 
between these 10 structures was 1.1 2 A. 



and oriented such that the H4a proton and the CH3 
group faced into the minor groove, whereas the H7a, 
H8a and H9a protons faced into the major groove 
(Figure 7). The IQ ring was angled by ~15° with respect 
to the modified guanine, but otherwise remained largely in 
plane with the damaged base. The helix was unwound 
between and X^, with a reduced helicoidal twist of 
30°. This was partially compensated by an increased 
twist of 9° between X^ and C . At base pair X^:C^^, the 
roll of the X^ purine decreased by 24°. This was 
compensated at base pair C**:G'^, where the roU decreased 
by 12°. Consequently, the A'^^-dG-IQ adduct induced a 
bend of 10° to the duplex. The IQ ring exhibited 
stacking with the flanking base pairs (Figure 8). IQ was 
stacked between G'^ and G" of the complementary 
strand of the C^:G'' and C*^:G'^ base pairs. The comple- 
mentary nucleotide, C'^, extruded into the major groove 
and did not exhibit stacking with the neighboring base 
pairs. The base opening between X^ and C'** increased 
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CORMA Value for Individual Bases 



0.8O 

oc" 0.SO 
0.4O 
0.3O 
0.20 
0.10 
0.00 



3 4 5 6 7 S 

Nucleotide 
CORMA Value for Individual Bases 



10 



12 



1.00 
0.90 
D.fiO 
D.70 
0.60 
CC 0.M 
0.4O 
0.30 
0.20 
0.1O 
D.OO 



24 



23 



22 



21 



20 



19 



18 



17 



Ifi 



15 



14 



Nucleotide 

Figure 6. Sixth root residuals (R^,) calculated using complete relaxation matrix calculations from the average of 10 structures emergent from the 
rMD calculations of the A^^-dG-IQ modified duplex. The black bars represent intra-nucleotide sixth root residuals, and the gray bars represent inter- 
nucleotide sixth root residuals. (A) Nucleotides C'-C'^ in the modified strand. (B) Nucleotides G'^-G^'' in the complementary strand. 



Table 4. Structural statistics for the N -dG-IQ modified duplex 



Average structure (obtained from 10 structures) 


RMS pairwise difference between structures 
RMS difference from average structure 


1.12 
0.75 




CORMA analysis for average structure" 


Intranucleotide 


Internucleotide 


Total 


R^^ 0.085 
Average error" 


0.083 


0.084 
0.0037 



"The mixing time was 250 ms. ^ , , , , , , , , , , , , -, 
"R'. IS the 6'^ root R factor: E[((Io);'^) - ((ic),''^)/ E((Io),"')]. 

''Average error: X!(Ic — lo)/", where the l^ are NOE intensities calculated from the refined structure, the lo are experimental 
NOE intensities. 



by 76°. This disrupted Watson-Crick hydrogen bonding. 
The other base pairs maintained Watson-Crick hydrogen 
bonding. 

DISCUSSION 

The A^^-dG-IQ DNA adduct has been of interest foUowing 
reports that it is more persistent than the C8-dG-IQ 
adduct in rodents and primates that were fed IQ in their 
diet (63). The synthesis of this adduct into oHgodeoxynu- 
cleotides (67) has allowed the conformation of the A^^-dG- 
IQ adduct at the G^ position of this sequence to be 



determined. This is a hot spot for two-base frameshift 
deletions in bacterial mutagenesis assays (68-71,73,74). 
In addition, human DNA polymerase (hpol) t\ produces 
two-base deletions when replicating past the A^^-dG-IQ 
adduct at the reiterated G^ position of the Narl 
sequence, in vitro (72). 

Conformation of the A^^-dG-IQ adduct 

The IQ ring intercalates when the A^^-dG-IQ adduct is 
positioned at the frameshift-prone G^ position of the 
Narl sequence (Figure 7). The strong NOE intensities of 
the IQ H4a and CH3 protons to the and C*^ HI' 
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Figure 7. Expanded view of the average structure calculated from 10 
structures emergent from the rMD calculations of the A'~-dG-IQ{X'') 
modified duplex, showing base pairs C<':G''', X'lC'** and C^.G". The 
view is from the major groove. The modified base pair X':C'* is shown 
in magenta, with the IQ moiety shown in blue. 



protons (Table 1) indicate that these protons face into the 
minor groove and establish the conformation about the 
bond between N^-dG and C5 of the IQ moiety. In 
contrast, NOEs involving the H8a proton of the IQ ring 
are primarily to bases G' , C'^ and G'^ of the complemen- 
tary strand (Table 1). The chemical shifts of the IQ H7a, 
H8a and H9a protons are observed between 6.5 and 
8.0 ppm, which is 1.3-2.0 ppm upheld as compared to 
the A'^^-IQ-dG nucleoside. This is consistent with the 
intercalated conformation and stacking of the IQ ring 
below the 5'-neighboring G" of the cornplementary 
strand and above the 3'-neighboring C^:G' base pair 
(Figure 8). Chemical shift perturbations corroborate the 
NOE data (Figure 4). The IQ H4a proton resonance, 
observed at 9.6 ppm, is 0.4 ppm upheld from the resonance 
observed for the modified A'^^-dG-lQ nucleoside (67), 
consistent with its location below G'^ and above C** 
(Figure 8). The IQ moiety displaces the complementary 
C'** base from the duplex, and flips it into the major 
groove. This is supported by smaller perturbations in 
chemical shifts for the H4a and CH3 protons as 
compared to the H7a, H8a and H9a aromatic protons of 
IQ. The displacement of the modified nucleotide 
toward the major groove (Figure 8) is supported by the 
downfield chemical shift change of 0.4 ppm for the H8 
and HI' protons of the modified base. The H6 proton 
resonance also experiences a downfield shift of 0.3 ppm. 
The stacking interactions of the IQ ring with the flanking 
bases C**, G ' and G" are reflected in the thermodynamic 
analysis of the adduct, in which the thermal melting tem- 
perature of 63° C is unchanged from that of the unmodi- 
fied duplex. 

Comparison to the A^^-acetylaminofluorene-dG adduct 

The other A^^-dG arylamiiie adduct that has been sub- 
jected to conformational analysis, although not in the 
Narl sequence of interest herein, is that arising from 
A^-acetylaniinofluorene (AAF; Chart 1) (76). The TV^-dG- 
AAF adduct conformation has also been examined using 
computational approaches (96). Zaliznyak et al. (76) have 
shown that the AAF moiety resides in the minor grove 
with its long axis directed toward the 5'-end of the 
modified strand. This shields the hydrophobic AAF ring 




Figure 8. Expanded views of the average structure calculated from 10 
structures emergent from the rMD calculations of the A'^-dG-IQ(X') 
modified duplex. Base stacking of the modified X':C'* base pair with 
the 5'-neighbor and 3'-neighbor base pairs. (A) Stacking of C^:G" 
above X':C"*. (B) Stacking of X':C"* above C^:G". Expanded views 
of the average structure of the corresponding C8-dG-IQ adduct at X'. 
(75) (C) Stacking of C'':G''' above X':C'*. (D) Stacking of X':C"* 
above C'*:G". In each instance, the modified base pair X':C"* is 
shown in magenta, with the IQ moiety shown in blue. 



from water. Similar to the -dG-IQ adduct, the modified 
nucleotide maintains the anti-conformation about the 
glycosyl bond. Notably, the A^^-dG-AAF adduct increases 
the stabihty of the DNA, which has been attributed to a 
favorable entropic effect (76). The present data reveal 
that the base-displaced intercalated conformation of the 
A^^-dG-lQ adduct at position G^ of the Narl sequence 
differs from that of the A^^-dG-AAF adduct, suggesting 
that the conformations of A^^-dG arylamine adducts vary 
rather than following a common motif. At the molecular 
level, the factors governing whether planar aromatic mol- 
ecules such as AAF or IQ favor DNA groove binding 
versus intercalation are not well estabhshed, but may be 
influenced both by their electronic structures and their 
respective geometries (97). Replication bypass studies 
have revealed that the A^^-dG-AAF adduct largely 
blocked DNA synthesis, but with some bypass and mis- 
incorporation of dATP opposite the lesion (98). 

Comparison to the C8-dG-IQ adduct 

When the C8-dG-IQ adduct was placed into the Nar\ 
sequence at the frameshift-prone G^^ position, the IQ 
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ring also intercalated into the duplex and the complemen- 
tary C'^ base was extruded into the major groove. The 
conformation of the C8-dG-lQ adduct also was 
characterized as base-displaced intercalated (75). Thus, 
at the G^^ position within the Narl sequence, both the 
C8-dG-lQ and A^^-dG-lQ adducts share a motif in which 
the IQ ring intercalates and C'^ is extruded into the major 
groove. However, the two conformations are distinctive. 
Apart from the difference in the regiochemistry of alkyl- 
ation (C8 versus A^^; Scheme 1), a major difference 
between the C8-dG-lQ and 7V^-dG-IQ adducts is that 
the C8-dG-IQ-niodified guanine adopts a syn conform- 
ation about the glycosyl bond, whereas the A^^-dG-IQ- 
niodified guanine maintains the anti-confonnation about 
the glycosyl bond (Figures 8 and 9). In addition, for the 
C8-dG-lQ adduct, rotation of the glycosyl bond into the 
syn conformation places the Watson— Crick hydrogen 
bonding edge of the modified dG into the major groove. 
The imino and amino protons are exposed to solvent. 
For the C8-dG-lQ adduct the orientation of the IQ ring 
with respect to the base is opposite to that of the A^^-IQ 
adduct, such that the IQ CH3 group and H4a and H5a 
protons face the major groove rather than the minor 
groove (75). The orientation of the C8-dG-IQ adduct in 
the duplex rotates the bulk of the IQ aromatic ring away 
from the flanking bases, resulting in a loss of base-stacking 
interactions, as shown in Figure 8. In comparison, the 
A^^-dG-IQ adduct appears to have more favorable 
stacking interactions with G'^. These differences may 
lead to differential processing during both DNA repair 
and DNA replication. 

Structure-activity relationships 

The A^^-dG-IQ adduct is less efficiently removed from 
genomic DNA by nucleotide excision repair (63,99). The 
NER machinery is thought to recognize bulky DNA 
damage that is destabilizing and distortive to the duplex 
(100-103). It has been proposed that the thermal stabil- 
ization of the A^^-dG-AAF adduct hinders NER (76). We 
observe that the of the A^^-dG-IQ adduct at position 
G^ within the Narl sequence does not destabilize the 
duplex (Table 1), correcting our original report (67). The 
Ta-i of the A^^-dG-IQ modified duplex is 63° C, and does not 
differ significantly from the unmodified duplex. This is 
remarkable given that the intercalated IQ moiety 
disrupts Watson-Crick hydrogen bonding and that the 
complementary C'^ base is displaced into the major 
groove. The stability of the iV^-dG-IQ modified duplex 
likely arises from favorable stacking between the IQ 
moiety and the neighboring base pairs (Figure 8). It is 
also interesting to note that unlike the A^^-dG-IQ 
adduct, the C8-dG-IQ adduct, which does not stack with 
the neighboring bases as well at this position (Figure 8) 
thermally destabilizes the duplex, reducing the by 4°C. 
Yeo et al. (104) examined AAF and AF C8-dG adducts 
within the Narl sequence and observed a correlation 
between the degree of destabilization induced by the 
lesions, binding affinities to the damage recognition 
protein XPC-RAD23B and overall NER efficiencies. 
Likewise, Zaliznyak et al. (76) attributed the increased 



stability of the A^^-dG-AAF adduct to its orientation 
within the minor groove and the entropy-favored release 
of waters from the duplex. Similar conclusions were 
reached by Cai et al. (105) who correlated thermodynamic 
stabilities and van der Waals interaction energies with 
repair efficiencies for stereoisomeric intercalated A'^-dA 
PAH adducts. Their studies showed that intercalated 
adducts with fewer DNA structural distortions and 
increased van der Waals interactions with neighboring 
bases correlated with reduced repair efficiencies. The 
HCA PhIP adduct has been compared with the cis- 
B[a]P-A^"-dG adduct in duplex DNA and in a nucleotide 
deletion duplex, and it was concluded that local stabiliza- 
tion induced by these adducts governs the ability of the 
(3-hairpins of NER proteins to recognize the damage 
(^106). In summary, it seems plausible that the thermal 
stability of the -dG-IQ adduct may, in part, explain 
the persistence of the A^'-dG-IQ adduct in rats and 
primates. 

If not repaired, the A^"-dG-IQ adduct is anticipated to 
be genotoxic. Indeed, IQ is an order of magnitude more 
mutagenic than is aflatoxin Bi in Ames assays. The muta- 
tions occur primarily at G:C base pairs (20,21). The rep- 
lication of the iV^-dG-IQ-adduct within the Narl sequence 
is influenced by the identity of the DNA polymerase. 
Because the damaged guanine remains in the anti-con- 
formation about the glycosyl bond (Figures 7 and 8), 
one might anticipate that the A^^-dG-IQ lesion should 
block Watson-Crick base pairing with incoming dNTPs 
during lesion bypass. Stover et al. (107) incorporated 
the 7V^-dG-IQ-adduct into the G' and G^^ positions of 
the Narl sequence and examined replication of the 
ohgodeoxynucleotides with Escherichia coli polymerases 
(pol) I (exonuclease deficient Klenow fragment), exonucle- 
ase deficient pol II and the Solfolobus solfatariciis P2 DNA 
polymerase IV (Dpo4), in vitro. At the G'' position, the A^"- 
dG-IQ adduct blocked the E. coli polymerases. Pol II exo~ 
favored correct incorporation of dCTP over dGTP but 
was unable to extend either of these initial insertion 
products. In contrast, the Dpo4 polymerase bypassed 
the A^^-dG-IQ adduct and produced an error-free 
product. The present studies do not necessarily predict 
the structure of the A^^-dG-IQ adduct during trans-lesion 
bypass. Consequently, it will be of interest to prepare 
complexes of bypass polymerases with A^^-dG-IQ 
modified template:priniers in an effort to determine how 
the A^^-dG-IQ adduct is accommodated during lesion 
bypass and how polymerases, e.g. the Dpo4 polymerase 
(108) allow bypass of this lesion. 

Bypass of the A'^^-dG-IQ adduct has been reported to be 
dependent upon its position in the Narl sequence. Choi 
et al. (72) have demonstrated that the human DNA poly- 
merase (hpol) r| can extend primers beyond template 
A^^-dG-IQ adducts. Pol r| correctly inserts dCTP and in- 
correctly inserts dATP. Analyses of hpol ri extension 
products reveal that a —2 bp deletion occurs with the G"^ 
A^^-dG-IQ adduct. In contrast, at the G' position replica- 
tion past the A^^-dG-IQ adduct results in error-free incorp- 
oration of dCTP, but further extension is inhibited and the 
polymerase stalls. In contrast, hpol ri does not yield —2 bp 
deletions with the C8-dG-IQ adduct located at position 



3460 Nucleic Acids Research, 2014, Vol. 42, No. 5 



G^. While further studies will be necessary to probe the 
basis for these observations, it is of interest to note that 
the stability of the A^^-dG-IQ adduct placed opposite a 
2-bp deletion increases as compared to the fully comple- 
mentary duplex, suggesting that the adduct may stabilize a 
2-bp strand slippage intermediate (67). At the G' position, 
the A^^-dG-IQ adduct is bypassed and extended by the 
E. coli polymerases and the Dpo4 polymerase, and 
error-free product is observed. Thus, it will also be of 
interest to examine the structure of the A^^-dG-IQ 
adduct when positioned at position G' of the Narl 
sequence. 

Summary 

Analysis of the A^^-dG-IQ adduct placed at position G"* of 
the Narl sequence (68-71,73,74), where it has been 
observed to cause —2 bp deletions when bypassed by 
hpol ri (72), reveals that it adopts a base-displaced 
intercalated conformation in which the H4a and CH3 
protons of the IQ ring face the minor groove and the 
H7a, H8a and H9a protons face the major groove. The 
IQ ring is shielded from water and stacks with the 5'- and 
3'-neighbor base pairs. Remarkably, despite this conform- 
ational perturbation, the A^^-dG-IQ adduct does not 
destabilize the duplex, which may correlate with the ob- 
servation that it is refractory to repair by NER (63,99). In 
addition, the IQ moiety disrupts the potential for Watson- 
Crick hydrogen bonding with incoming dNTPs, which 
perhaps explains why this lesion blocks DNA synthesis 
by many polymerases. 
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